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ABSTRACT 



The effects of heat transfer on the kinetic energy tendencies of a 
perturbation superimposed on a basic zonal current are investigated. The 
rate at which kinetic energy is gained in both the adiabatic and diabatic 
models increases with time for medium wave lengths, but the rate of in- 
crease is less for the diabatic case. At the longer wave lengths the 
energy change for the diabatic model shows a decrease with time becoming 
negative toward the end of the time interval investigated. The adiabatic 
model, in contrast, has zero energy tendency at long wave lengths. The 
rate of kinetic energy change increases as larger values of thermal wind 
are selected, and for the adiabatic model the wave length range for which 
energy changes occur also expands with increasing thermal wind. The 
energy gain initially is a maximum when the thermal wave lags the pressure 
wave by a phase angle of 90 degrees and the maximum tendency to lose 
energy occurs when the lag is 270 degrees. 



UISfJlRl 

1IAVAL POSTGRADUATE SCHOOL 
LICHTETEY, C'LTF. 03010 

TABLE OF CONTENTS 

Section Page 

I. Introduction 7 

II. Development 9 

III. Discussion of Results 16 

IV. Acknowledgements 22 

V. Bibliography 23 



3 



LIST OF ILLUSTRATIONS 



Figure 




Page 


1 . 


Amplifying (+D) or dampening (-D) factor as a function 
of wave length for thermal wind lOm/sec - No Heating 


24 


2. 


Amplifying (+D) or dampening (-D) factor as a function 
of wave length for thermal wind 20m/sec - No Heating 


25 


3. 


Amplifying or dampening characteristics B+D, B-D, B, 
+D, -D for thermal wind lOm/sec - Heating 


26 


4. 


Amplifying or dampening characteristics B+D, B-D, B, 
+D , -D for thermal wind 20m/sec - Heating 


27 


5. 


Wave Speed - Heating (A H) vs No-Heating (A NH) for 
thermal wind 20m/sec 


28 


6. 


KE tendency (ergs/sec) for phase lag 60,90,150,240, 

270 and 330 degrees for thermal wind 15m/sec - No Heating 


29 


7. 


KE tendency (ergs/sec) for phase lag 60,90,150,240, 

270 and 330 degrees for thermal wind 15m/sec - Heating 


30 


8. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 0 hours and phase angle 90 degrees - No Heating 


31 


9. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 0 hours and phase angle 90 degrees - Heating 


32 


10. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 6 hours and phase angle 90 degrees - No Heating 


33 


11. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 6 hours and phase angle 90 degrees - No Heating 


34 


12. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 12 hours and phase angle 90 degrees - No Heating 


35 


13. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 12 hours and phase angle 90 degrees - No Heating 


36 


14. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 24 hours and phase angle 90 degrees - No Heating 


37 


15. 


KE tendency (ergs/sec) for thermal winds 10,15 and 20 
m/sec at 24 hours and phase angle 90 degrees - Heating 


38 



4 



TABLE OF SYMBOLS 



J = 

V = 

io = 

V = 

3 - 

i = 

T = 

to = 

s = 



fi - 
a - 

j = 




r - 

a = 

L- 



t = 

t = 

u = 

%*= 



the relative vorticity 
the wind velocity vector 

the average value of the Coriolis parameter and equal to the 
the value at latitude 45N 

the gradient operator in pressure coordinates 

the acceleration of gravity 

the height of the pressure surface 

the temperature 

the total derivative of pressure with respect to time (vertical 
velocity) 

/ 3T_ 

the static stability parameter p Op ) 
the Rossby parameter 
the gas constant for dry air 
the Jacobian operator 




a measure of static stability -p S R. 
the specific heat for constant pressure 



the stream function ( 2.^ 

the thickness function ( -3- 
as 7T = 3.14159 V l* 

the heat gain or loss 

the atmospheric pressure 
/ 

the wave number f J 

the wave length 

the phase speed of the waves 

time 

the basic zonal wind 
the basic thermal wind 




until specified 
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I. Introduction 



To date simple baroclinic models have been unable to improve signif- 
icantly on the results obtained by one parameter models in the prognosis 
of mid-tropospheric circulation- Recent investigations have shown that 
the incorporation of previously neglected terms, which were felt to be 
negligible for short range forecasting, can significantly alter the 
characteristics of synoptic waves produced by baroclinic forecasting 
models, thus bringing these waves into closer proximity with snyoptic 
waves observed in nature. 

4 In one of these investigations, Haltiner to considered the influ- 
ence of sensible heat exchange on the properties of harmonic perturba- 
tions superimposed on a zonal current- The study revealed that this 
type of heating reduces the instability of short and medium waves, but 
that the instability of long waves ia increased. 

The purpose of this work is to compare the kinetic energy changes 
averaged over the interval of one wave length in a baroclinic model with 
heat transfer present to one when heating is absent. In particular, the 
manner in which heating alters the instability of long waves will be 
examined. A simple two-level quasi-geostrophic model is constructed * 
and an equation expressing the time rate of change of kinetic energy 
in a harmonic perturbation superimposed on a zonal current is derived- 

9 

The kinetic energy is computed numerically for the range of wave lengths 
corresponding to sy.noptically observed waves. At each wave length com- 
putations are made for varying values of thermal wind and phase lag be- 
tween the pressure wave and the temperature wave. The change of kinetic 
energy with time is evaluated in time steps of three hours over a thirty 
hour period and the value of the tendency for the heating and the non- 
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heating model are compared. 

In addition, computations are made of the values of C^/L, which is 
proportional to the growth rate of the wave, and the values of phase 
speed to verify that the results obtained from the simple model used are 
in qualitative agreement with those obtained by Haltiner. 

Interest throughout the study is confined to unstable waves only 
and solutions which lead to neutral waves are eliminated. 
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II . Mathematical Development 



The atmosphere is represented by a four-layer baroclinic model 
using the 750 and 250 millibar levels as data levels and centering atten- 
tion at 500 millibars. Since the effects of latitude variation are well 
known, latitude is held fixed at 45 N. The vertical velocity is assumed 
to be zero at both the top and the bottom of the model atmosphere, the 
effects of friction at the earth's surface are neglected, and the static 
stability of the atmosphere is taken to be, constant. 

The vorticity equation is applied at the levels corresponding to 
250 and 750 millibars in the following form: 






+- 



St 



v-vOjO= ]° v-v. 



( 1 ) 



The continuity equation is used to substitute for the divergence in the 
above equation and then the derivative of the vertical velocity with 
respect to pressure is approximated by a finite difference ratio. Geo - 
strophic approximations are substituted for the vorticity and wind. The 
sum and difference of this equation at the data levels are computed with 



the following results 

X 



V 



(2 



o 



At the middle level (500 mb) the First Law of Thermodynamics is 
applied in the following form: 

- l * V • V T * W 

3 1 r 

where / cv \ 

\ “ \ 3 '\> c f X ♦ 
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The temperature at this level is approximated using the equation of 



state and the hydrostatic equation to give 




This approximation is used in (4) and used to replace Co in (3). 



With the following substitutions: 




our system of equations becomes 

Wsi) *J(^V^) + jj(Tr J VV) + f 

- T K^4>) 




(5) 



( 6 ) 



Let the basic flow and the surface be considered to be in thermal 

equilibrium. Similar to Doos dl , the exchange of sensible heat is 

then assumed to be proportional to the difference between the surface 

perturbation air temperature (T ) and the temperature of the ground (T ) 

® £ 

and a decreasing power of pressure in the following manner: 

f - 

where G is a proportionality constant. 

Since ridges in the thermal wave will correspond to heat sinks for 
the air, and troughs will be heat sources, let the transfer of heat be 
proportional and opposite in sign to the perturbation thickness. There- 
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fore 



KTT 



/ 



A. ^ 
"‘VI* s 



( 8 ) 



where K is a proportionality constant with units sec“^. 

To linearize the system (5) and (6), assume that the solutions are 
of the form 







and that 






Dropping terms involving products of perturbation quantities leads 



to a homogeneous linear system MN=0, where M is the matrix 

-a* 2 * 

f 



iklfi* ^ i 

LkiL*y-**) 



bk, (t-V)(u\k l )+Lkj2yU V Kj 



(9) 



and N is a column matrix of ifo and . In order to have non-trivial 
solutions to this system of equations the determinant of the coefficient 
matrix M must vanish. Setting this determinant equal to zero yields the 
frequency equation 




4 ) 



where 




It can be seen from this equation that the presence of heating in- 
sures that the phase velocity will be a complex number regardless of 
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the wave length. If the heating is taken to be zero, it can be shown 
that the phase velocity will be complex only for a certain range of wave 
lengths. The width of this range will depend on the strength of the 
thermal wind and the value of the static stability. 

Using the primary roots of the complex number under the radical sign, 
the phase velocity has two solutions: C + when using the positive sign in 
front of the radical, and C_ when using the negative sign. The complete 
solutions for y and 7T are therefore given by 




ik(x-i+t) 

JL 



-h 



, - ik U-*-*) 

Vo *- 



( 11 ) 



and 

T- K *- +■ T 0 -*• 



(12) 



where the primes are omitted and henceforth y and TT will be perturba- 
tion quantities unless otherwise specified. 

Next take as initial conditions at time t=0 




7 r - B*. 



where A and B are real numbers and ^ is a relative phase difference 
which is positive when the thermal field lags the stream field. From 
(11) and (12) it follows immediately that at t=0 



X + TT = (13) 

Along with (13), use of the linearized form of (5), once for C + and 
once for C_, yields four equations in ^ , (j/ t , 77^ , and 77^ . Letting 
Cp equal the Rosshy wave speed, the following solutions result: 



12 



(14a) 



j4 r * A - 



[* fa-**)- f] 



K'O 



p±^±_e£ | 7 

^ <V‘J 

f (c.-O 



(14b) 



__f -o t( * 

T7^ - 8-e. - 



(14c) 



7 r 



(14d) 



r(t*-o „ 



The above solutions for r- ^ and are complex func- 

tions. Therefore, they can be separated into their real and imaginary 
parts 



K - C * * 

and similarly for and 77^ . 



As mentioned previously, the phase velocity can be a complex num- 
ber depending on whether heating is present and the values of the thermal 
wind and static stability. The phase velocity can then be represented by 






C l_ =■ a + c b - 4-id) 



(15) 



where 



- / /,£ . / "j 



a ■ - 



and 






a.iJ : 1)7. *fUi* ) 

* (ft. **£.*+**)* J f y 



'/z. 



(16) 
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The kinetic energy of the perturbation may be expressed in terms of 



the square of the velocity of the perturbation wind. Since p is inde- 
pendent of y and considering a unit mass, the kinetic energy per unit 
mass is 






2 . 



_ l 

" L 




A solution for V is obtained by substituting into (11) the solutions 
for C + and C_ as given in (15) and then differentiating (11) with respect 
to X. The real part of V is squared, multiplied by one half (1/2) and 
integrated over an interval of one wave length to give the average 
kinetic energy of the perturbation over one wave length as follows: 

ijVk . ^ +;;) , 

* -n- k [/%& fc - c t'z A- ay)] 

4- rrk ^ hr ] 



where Y = and 7T in this and subsequent equations is 7 7" = 3.14159. 

The time rate of change of this kinetic energy is then obtained by 
differentiating (17) with respect to time with the following result :: 



__TL*7 (b ♦ Vr )] 

+ C) ir)] 



(18) 



Y -e- lbH j il(b O*^)] 
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+ ^V bk< 



[u;£:0> 



am 2\ f-l-G CM 



Jr) - i is ft (t> 



C*}}#- 4 z &m 



»)] 



This energy change equation when evaluated at some time t=t^ repre- 
sents the tendency of the perturbation to gain or lose kinetic energy. 
This gain or loss then shows whether the wave is amplifying or dampening. 
Since only unstable waves are studied, a solution leading to zero corres- 
ponds to a neutral wave. 

Haltiner [3] has shown that the mean energy changes over one wave- 
length for adiabatic conditions can be expressed as 



where KE and PE are the mean kinetic and potential energies. The quantity 



into (12) and solving for Ca> from (4), the energy transformation func- 
tions on the right hand side (19) may be determined. The kinetic energy 
changes that have been examined in this paper are thus a consequence of 
this energy transformation from potential to kinetic energy. 




and 



(19) 





equal to 7f ; so after substituting the values for C + and C 
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III. Results 



It is obvious from equation (18) that a determination of precisely 
how heating affects the stability characteristics of harmonic waves would 
be difficult. Examination of the frequency equation (10) does show that 
with heating present an instability term (term b) is always present, and 
also that heating leads to an instability term evolving from the radical 
part of the equation (term d) . Term b always leads to a dampening effect 
while term d may have either dampening or amplifying characteristics. 
Therefore, the combined effects of these two terms will determine the 
energy change characteristics of the heating model. In the adiabatic 
model the energy change is dependent on only term d, which takes on a 
different value in general when heating is present 

To arrive at quantitative results for energy changes in the two sit- 
uations, representative values of the different parameters involved are 
selected and equation (18) is solved numerically. The latitude is fixed 
at 45 degrees north latitude, the basic flow zonal wind remains constant 
at 20 m/sec, and the static stability also remains constant with = 

-6K/100mb and sojt^ = 2.4 x 10 -1 ^ cnT^. The initial amplitude of the 

stream function and the thickness function are taken to be equal with a 
10 2 

value of 20 x 10 cm /sec which corresponds to a z-amplitude of 20km. 
Since both the thermal wind and phase lag are of significant importance 
in the energy change, four values of thermal wind (5, 10, 15, and 20 
m/sec.) and 12 phase lags, 30 degrees apart, are used. 

It remains only to select a representative value for the heating 
proportionality constant. The value for G in (7) is chosen to be 

IOO CR.6,% °K~' <,ec~‘ 
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and the value for r is r=2. This is similar to values used by Doos 



» • 



D] 



and corresponds to a total heat transfer of 63 cal per square centimeter 
per day to the entire column of air for every one degree temperature dif- 
ference between the ground and the air at the surface. To determine a 
value for K, the initial perturbation thickness is translated into a 
temperature perturbation by use of the hypsometric equation. Since the 
surface air and the ground are assumed to be initially in thermal equil- 
ibrium, this temperature perturbation also equals the maximum difference 
between the air and the ground. A perturbation thickness of 20 x 10 cm / 
sec which is the selected initial value translates into a temperature of 
6.54K and corresponds to a heat transfer of 654 ergs g~^ sec“^. A value 
for K can now be calculated from (8). The value for K at 500 millibars 
is then arrived at by using r=2 and turns out to be 



Using the above values, the first items determined were the stabil- 
ity characteristics of the waves. As mentioned in the introduction this 
is the value of the complex part of the phase velocity multiplied by k, 



1 is for no heating and a thermal wind of lOm/sec. Figure 2 is again for 
no heating and now the thermal wind is 20m/sec. Figures 3 and 4 are the 
corresponding situations for the heating model. 

To gain sufficient resolution between the values for each model at 
interested wave lengths the graphs were scaled such that the values at 
short wave lengths often exceeded the bounds of the graph. Where the 
bounds have been exceeded such that the graph plot is a straight line 
parallel to the x-axis, the actual values at these wave lengths are 
greater than those depicted on the graph. In cases where the values 
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are significant to the problem, special mention will be made of the values 
It should also be mentioned that the data points are obtained at intervals 
of 1000 kilometers and between data points the computer simply makes 
straight line extrapolations. 

Examination of figures 1 through 4 reveals the following items. Heat 
ing leads to a smaller amplification factor than no heating and a larger 
dampening factor. The values at 5000km, which are off the graph for a 
thermal wind of 20m/sec, also show that the no heating leads to a larger 
amplification and less dampening. The effect of thermal wind is readily 
seen to increase the instability characteristics even though the value of 
the b term in the heating model is independent of thermal wind. Also it 
is important to note that in the adiabatic case the instability is con- 
fined to a finite wave length interval which depends on the thermal wind 
and static stability. In contrast the presence of heating leads to an 
instability factor at all wave lengths, however the amplification ap- 
proaches zero with increasing wave length. These results are in agree- 
ment with similar computations made by Haltiner. 

Figure 5 depicts the relationship between the speed of wave propaga- 
tion in the two models. Ignoring the value at wave length zero for the 
heating case, heating leads to a decrease of phase velocity of the wave 
at all wave lengths. This difference is, however, quite small at 4-6000 
km. This also is in general agreement with results obtained by Haltiner. 

Computations of the time rate of change of kinetic energy show some 
interesting features. First of all, if the energy change tendency at 
time t=0 is evaluated for several different values of phase lag, the max- 
imum positive effect occurs at a phase lag of 90 degrees and the maximum 
negative change is at a phase lag of 270 degrees. The pattern is sym- 
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metrical about these two points for both models with the exception that 
for the adiabatic case the values are identically zero at phase lags of 
0 and 180 degrees while in the case of heating they are small but not 
identically zero. A plot of the kinetic energy tendency at time t=0 for 
several different phase angles is shown in figures 6 and 7. As the time 
was increased from zero the phase lag of maximum positive and maximum 
negative energy change tended to increase, so that by 12 hours several 
situations (wave length and thermal wind values) showed a maximum posi- 
tive value at 120 degrees of lag and maximum negative at 300 degrees. 
Also, in all situations the magnitude of the positive rate of change was 
greater than that of the negative change. Again the graphs indicate that 
the presence of heating leads to kinetic energy changes at all wave 
lengths while its absence restricts the energy change to a finite inter- 
val. This phenomena remains true throughout all situations investigated. 

Another point that should be brought out at this time is that in 
the derivation of equation (18), both wave components of (11) and (12) 
were used. The energy change equation, therefore, represents the com- 
bination of the damped and amplifying components. Hence, the values of 
wave length, thermal wind, and phase lag which result in positive energy 
change represent situations when amplification will take place and vice- 
versa. 

When computations for the kinetic energy tendency for the heating 
and non-heating case were compared at t=0, the values were identical 
for corresponding situations (with the exception of phase lag 0 and 180 
degrees which was previously mentioned). This seemed quite surprising 
at first since equation (17) for the diabatic and adiabatic cases is 
quite different. Upon further reflection it is apparent that since 
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heating has had no time to exert an influence, the values would be the 
same. Figures 8 and 9 show the tendency at t=0 for several thermal winds. 
On the curves for a thermal wind of 20m/sec, the values at wave lengths 
5000, 6000 and 7000 kilometers can be seen to be equal in both cases. 

Next the tendency is evaluated at t=6 hours. In both cases the 
values increased in magnitude but now an examination of figures 10 and 
11 show that the values for the heating case are smaller than those for 
no heating. When time is allowed to increase further, the difference in 
tendency values for the two cases becomes larger. Figures 12 through 15 
are for t=12 and t = 24 hours and show how the difference increases. By 
24 hours the values at 5000 kilometers are off the graphs but the raw 
data still gives a higher value for the adiabatic situation. 

Values of energy change were actually tabulated out to 33 hours at 
a phase lag of 90 degrees. In the adiabatic model the tendency value 
increased out to 33 hours in all cases. When heating was present the 
value also increased with time for medium wave lengths with the exception 
that after 30 hours for a wave length of 5000 kilometers and a thermal 
wind of 5 and lOm/sec, the value started to decrease. This suggests that 
not only does the tendency increase at a slower rate when heating is in- 
troduced, but also that it reaches its maximum value sooner. At the 
longer wave lengths the value of the tendency decreases with time such 
that by 24 hours kinetic energy is actually being lost by the waves with 
wave lengths larger than 10,000 to 11,000 kilometers. It was mentioned 
earlier that the phase lag of maximum energy change increases with time 
also. This change of phase lag does not, however, alter the results that 
heating leads to smaller values of tendency as the time increases. 

So far we have examined primarily the spectrum of medium length waves 
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to determine the effects of heat transfer on the energy changes. At the 
long wave lengths the addition of heat transfer to the model leads to 
only a small rate of energy change. This, however, is significant since 
in the adiabatic case the energy change is identically zero at long wave 
lengths. That heating insured the presence of an instability term at 
all wave lengths was observed in equation (10) and again in the results 
of figures 3 and 4. These did not, however, indicate whether the effect 
was to amplify or dampen the wave. Examination of the kinetic energy 
tendency due to heating clearly shows that the effect depends on the time 
and the phase lag. For example, the series of graphs on energy change 
show that out to at least 24 hours, if the phase lag is 90 degrees, the 
heating leads to an amplifying wave at long wave lengths. Beyond 30 
hours the heating leads to a dampening effect. 

The energy tendency at the very short wave lengths was not examined 
comprehensively and no explanation is offered for the very large values 
of energy change observed. It can be seen from the graphs, however, that 
besides displaying very large values, the kinetic energy tendency may 
change abruptly from a negative to a positive value or vice versa as the 
time or wave length varies. 
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